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ABSTRACT 

Crucial to a quantitative understanding of galaxy evolution are the properties of the inter- 
stellar medium that regulate galactic-scale star formation activity. We present here the results of 
a suite of star formation models applied to the nearby blue compact dwarf galaxies NGC 2915 and 
NGC 1705. Each of these galaxies has a stellar disk embedded in a much larger, essentially star- 
less H I disk. These atypical stellar morphologies allow for rigorous tests of star formation models 
that examine the effects on star formation of the H I, stellar and dark matter mass components, as 
well as the kinematics of the gaseous and stellar disks. We use far ultra-violet and 24 jiia imaging 
from the Galaxy Evolution Explorer and the Spitzer Infrared Nearby Galaxies Survey respectively 
to map the spatial distribution of the total star formation rate surface density within each galaxy. 
New high-resolution H I line observations obtained with the Australia Telescope Compact Array 
are used to study the distribution and dynamics of each galaxy's neutral inter-stellar medium. 
The standard Toomre Q parameter is unable to distinguish between active and non-active star 
forming regions, predicting the Hi disks of the dwarfs to be sub-critical. Two-fluid instability 
models incorporating the stellar and dark matter components of each galaxy, in addition to the 
gaseous component, yield portions of the inner disk unstable. Finally, a formalisation in which 
the H I kinematics are characterised by the rotational shear of the gas produces models that 
very accurately match the observations. This suggests the time available for perturbations to 
collapse in the presence of rotational shear to be an important factor governing galactic-scale 
star formation. 

Subject headings: galaxies: star formation — galaxies: ISM 



1. Introduction 

Galaxies are gravitationally bound systems of 
stars, gas, and dark matter. The interaction 
and evolution of these mass components results 
in a change of the observable properties of the 
galaxy as a whole. An understanding of the pro- 
cesses governing galaxy formation and evolution 
is still a major challenge to astronomers. The 
goal is to quantitatively define the relationships 
between galactic-scale star formation activity and 
the properties of the in ter-stellar medium (ISM) 
fe.g. IWvder et al.ll2QQ9[ ) . Current observing facili- 



ties allow for high-quality multi- wavelength obser- 
vations of comparable resolution and sens itivity to 



be ca rried out for nearby galaxies (e.g. Lerov et al 



20081 ) . These observations are used to test current 



theories of galactic-scale star formation. 

The cold gas component of a galaxy is its fuel 
for star formation. On global as well as localised 
length scales the distribution of these mass compo- 
nents is known to correlate with the star formation 
activity _('Kennicutt""l983l, Il989l Il998l : iBigiel et al.l 
20081 : Eeroy et al. 2008). In addition to the dis- 



tribution of the gas, it is also the kinematics that 
play a role in regulating the star formation activ- 
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ity. For example, while self-gravity drives gas per- 
turbations to collapse under mutual gravity, the 
centrifugal forces associated with their rotational 
motions will impede the structure growth by coun- 
teracting the in ward gravitational force (jToomre 



19641 iKennicutt 1989). Similarly, rotational shear 



will tear adjacent gas clouds apart before they 
have the cha nce to coalesce to fo rm a single denser 
body of gas ( Wvder et al.ll2QQ9[ ). The self-gravity 
of a galaxy's ISM is determined by the gravita- 
tional potential of its mass components. The stel- 
lar potential, for example, contributes to the self- 
gravity of the gas, thereby making it more suscep- 
tible to gravitationa l collapse and hence st ar for- 
mation (e.g. Rafikov 200llLer^elaD[2Q^ Star 
formation models are used to understand how all 
of these processes work together to control the star 
formation in galaxies. Both the ability and inabil- 
ity of these models to accurately describe the star- 
formation contributes to a further understanding 
of the complex interplay between the various ISM 
properties that regulate the star formation activ- 
ity. 

In order to compare modelling results to obser- 
vations, a quantitative measure of a galaxy's star 
formation activity is required. One of the most 
widely used tracers of high-mass star formation is 
the Ha emission line arising from the recombina- 
tion of ionised hydrogen. This traces a galaxy's 
star formation over the relatively short lifetimes 
(a few million years) of the most massive stars. 
In recent years, GALEX observations of nearby 
galaxies have allowed the ultra-violet (UV) flux 
originating from the photospheres of high-mass O- 
and B-type stars to be used as a diagnostic of the 
directly o b serva ble star formation rate (SFR). As 
iLee et aP (|2QQ9[ ) point out, a galaxy's UV flux 
probes a fuller mass spectrum of massive stars, 
and thus measures star formation averaged over 
a longer ~ 10^ yr time-scale. The UV emission 
alone is, however, an insufficient tracer of the SFR. 
Several authors (e.g. ISalim et al. 2007 and refer- 
ences therein) have compared Ha and UV SFRs 
for nearby galaxies. The general consensus is that 
UV SFRs tend to be slightly lower than Ha SFRs. 
The discrepancy is attributed to the effects of in- 
ternal dust extinction which lowers the observed 
UV flux. When UV emission from high-mass stars 
is absorbed by dust, it is r e-emit ted at infrared 
wavelengths. As iKennicuttI (| 19981 ) points out, in- 



ternal dust extinction serves as one of the largest 
sources of systematic error in SFR measurements. 
Combining a measure of a star-forming galaxy's 
infrared emission with a measure of its UV emis- 
sion yields a measure of the total SFR that is 
consistent with the rate inferred from Ha spectro- 
scopic imaging. fCalzetti et all (2007t ). for a sam- 
ple of 33 nearby galaxies observed with the Spitzer 
MIPS instrument at 24 /im, showed that using the 
24 /im emission to account for the dust extinction 
allows for an accurate determination of the total 
number of ionising photons in an HII region. 

Dwarf galaxies are particularly useful labora- 
tories in which to study star formation. Being 
morphologically and dynamically simpler systems 
than larger spiral galaxies, they allow us to strin- 
gently test star formation theories in the domains 
of low gas surface densities, extreme dark matter 
content, and lower rotational shear rates. This 
paper presents the results of a detailed study of 
the star formation activity in two nearby dwarfs, 
NGC 2915 and NGC 1705. Each galaxy has a 
small stellar disk embedded in a much larger Hi 
disk. Using the Australia Telescope Compact Ar- 
ray we have obtained new deep, high-resolution 
Hi line observations of these galaxies in order to 
study their extended Hi disks in unprecedented 
detail. We aim to link the Hi characteristics of 
each galaxy to its observed star formation activ- 
ity. Far ultra-violet and 24 /im imaging from the 
GALEX and Spitzer satellites, respectively is used 
to quantify the star formation activity. We com- 
bine all of these data to produce a suite of detailed 
star formation models that are used to understand 
the observed star formation in the context of the 
distribution and kinematics of the neutral ISM. 
T he results of our an alyses are compared to those 
of^ Lerov et al. ( 2008[ ) who carried out similar anal- 
yses for their sample 23 nearby late-type galaxies 
from The Hi Nea rby Galaxy Survey (THINGS, 
Walter et al1l2QQ8[ ). This is done so that the re- 



sults can be interpreted in the context of other 
more typical star- forming galaxies. 

The structure of this paper is as follows: Sec- 
tion [2] introduces NGC 2915 and NGC 1705. The 
new Hi, ultra-violet and 24 /im data sets that 
are utilised in this work are presented in Sec [3l 
The method of estimating the total star forma- 
tion rates is presented in Sec. |4] together with to- 
tal star formation rate maps of the galaxies. Each 



2 



of the star formation models is introduced and de- 
scribed in detail in Section[5l The modeling results 
for each galaxy are also shown in Section [H with 
comparisons to the results of iLerov et alJ (|2QQ8l ) 
included. The total star formation activity of each 
galaxy is used to estimate the possible unseen H2 
content of each system in Section [6l The results 
of this paper are summarised in Section [71 

2. NGC 2915 and NGC 1705 



At a distance of 4.1 ± 0.3 Mpc (jMeurer et al 



[2003) NGC 2915 is a nearby galaxy with the 
optical properties of a blue compact dwarf and 
the Hi characteristics of a late type spiral. The 
optical disk has a 5-band absolute magnitude 
of Mb = —15.90 and 25th magnitude isophotal 
radii of 1.2 kpc an d 1.9 kpc in the B - and R- 
bands respectively (Meurer et alJ Il994h , for the 
here adopted distance of 4.1 Mpc. The stellar 
disk is contained within a very small region at 
the centre of a large H I disk which extends out 
to > 22 ^-band s cale lengths (|Meurer e t al. 1996; 
Elson et al.ll2010[ ). The Hi disk has well-defined 
spiral structure as well as a bar-like feature co- 
located with the stellar disk. No significant star 
formation is observed in the outer Hi disk, only 
a few faint H ll regions have be en detected from 
deep Ra imaging (Meurer et al.iri999: .Werk et al.l 
2010h . 
Both 



Meurer et al. 



Elson et al. 



(|1996[ ) and 

(|201 0l) have used the extended H I disk of NGC 2915 
to trace its gravitational potential out to large 
galactocentric radii. Both investigations find the 
galaxy to be extremely dark-matter-dominated, 
with a total mass-to-blue-lig: ht ratio as hi^h a s 



Mtot/Ls ^ 140 Mq/Lb,q (|Elson et al.l l2010[ ). 



The Hi velocity field o f NGC 2915 was anal- 
ysed in further detail by Hi son et al.l (|201lf ) who 
find dynamical evidence for elliptical streaming 
and radial fiow patterns within the Hi disk, as 
well as a po s sibly aspherical dark matter halo. 



Bureau et al.l (|l999f ) showed that the Hi spiral 
structure can be caused by gravitational torques 
associated with a slow-rotating tri-axial dark mat- 
ter halo. 

NGC 1705, another blue compact dwarf, is 
well-known for hosting one of the most intense 
starbur sts (relative to its mass) in the local uni- 
verse. iTosi et all (|200lh used the HST to re- 



solve the brightest red giant stars in the galaxy, 
thereby estimating a distance of 5.1 ± 0.6 Mpc. 
The system has a 5-band absolute magnitude 
of Mb = -15.6 ± 0.2 (Marlo we et al.l Il999l ). 



Meurer et al.l (|l992l ) identified two stellar popu- 
lations in NGC 1705 which they refer to as the 
high and low surface brightness populations. The 
young stellar disk has a 5-band 25th magnitude 
isophotal radius of 1.2 kpc for the here adopted 
distance of 5.1 Mpc. The galaxy's intense star 
formation activity is concentrated in its central 
high surface brightness stellar population, and 
is driven mainly b y a powerful super star clus- 
ter, N GC 1705 - 1 (jSandagd [19781 : iMelnick et al 



1985f ). which contributes almost ha lf of the total 



ultra - vi olet luminosity o f the g alaxy (jMeurer et al 



19981 ). lAnnibali et al.l (|2003f ) found evidence for 
star formation activity that commenced at least 
5 Gyr ago. The authors also showed the star for- 
mation history of NGC 1705 to be very complex, 
confirming the existence of a recent burst of star 
formation between 10 and 15 Myr ago. Consis- 
te nt with the HS T observations of NGC 1705 - 
1, lAnnibali et al.l (j2003) found the young and 



intermediate-aged stars to be concentrated near 
the very centre of the galaxy. 

Like NGC 2915, the stellar disk of NGC 1705 
is embedded in a large, extended Hi disk. The 
first H I synthesi s observatio ns of the galaxy were 
carried out by M eurer et al.l (1998) who provided 
dynamical evidence for a galactic blowout that is 
powered by the central super star cluster. No dis- 
tinct spiral structure is visible in the H I distribu- 
tion. The general Hi kinematics of NGC 1705 
are those of a rotating disk. The dark matter 
properties of the g alaxy were first modelled by 



Meurer et al.l () 19981 ) who showed the system to be 



dark-matter-dominated at nearly all radii. 
3. Data 

Both galaxies have been observed at 21 cm, 
infrared, optical, and ultra-violet wavelengths. 
This section presents and describes the multi- 
wavelength data sets utilised in this paper. 

3.1. Hi line observations 



Meurer et al 



(jl996f ) and iMeurer et al.l (pi 
used the Australia Telescope Compact Array 
(ATC A) to carry out the first H I synthesis obser- 
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vations of NGC 2915 and NGC 1705, respectively. 
Meurer et al.l ( 1996[ ) used their data to produce 
the first mass models of a blue compact dwarf. 
In this work we utilise new H I synthesis observa- 
tions from the ATCA with significantly improved 
sensitivity and spatial resolution. 

The Hi li ne da ta for NGC 2915 are those of 
These data have a spatial 



Elson et al 



resolution of 17.0" x 1S.2'\ a fa c tor o f - 2 im- 



(|l996f ) observa- 



provement on the iMeurer et al 
tions. The r.m.s. of the noise in a channel is 
0.6 mJy beam~^, the chan nel width is 3 . 2 km s" 

(|2010h 



The reader is referred to Elson et al 



for 

a full description of the observations setups, the 
data reduction procedures, and a presentation of 
the various H I data products including the chan- 
nel maps. An H I surface densit y map of NGC 2915 



\i K 

produced using the data from l Elson et al.l (|2010l ) 



is shown in the top panel of Fig.[Tl For reference, a 
3.6 jam IRAC Spitzer image of the evolved stellar 
disk of the galaxy is shown in the bottom panel. 
The Hi observations are able to resolve the cen- 
tral region of the H I disk into two over-densities. 
Also visible is the clear spiral structure as well as 
a plume-like feature extending from the centre of 
the galaxy towards the north-west. It is this Hi 
surface density map that is used to investigate the 
star formation laws in this paper. 

The new NGC 1705 ATCA Hi synthesis ob- 
servations used in this work will be discussed in 
detail in our upcoming publication (Elson et al 
2012, in prep.). For this work it suffices to say 
that approximately 86 hours worth of on source 
observations were obtained for the galaxy using 
baselines up to 6 km in length to produce an Hi 
data cube with a spatial resolution of 16. 7^' x 14. 5^' 
and a channel width of 3.48 km s~^. The r.m.s. 
ffux of the noise in a line-free channel of the cube 
is ^ 0.7 mJy beam~^. An Hi surface density 
map produced from these data is shown in the 
top panel of Fig. [2l A 3.6 jam IRAC Spitzer im- 
age of the stellar disk is again shown in the bot- 
tom panel of the figure for comparison. The new 
H I data clearly resolve the central H I concentra- 
tion of the galaxy into two over-densities with a 
combined mass of ~ 3 x 10^ M©. These central 
Hi over-densities have their peaks separated by 
~ 0.8 kpc, and straddle the extremely luminous 
super star-cluster. 
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Fig. 1.— Top panel: NGC 2915 Hi sur- 
face density map produced using the data of 
Elson et al.i(|2010h . The colour scale is described 
in units of Mq pc~^ by the colour bar. The 
black ellipse represents the i?-band R25 radius 
of 1.9 kpc. The hatched circle in the lower left 
corner of the panel represents the half power 
beam width of the synthesised beam. Bottom 
panel: 3.6 jam IRAC Spitzer image of the old 
stellar disk of NGC 2915. The black ellipse is 
the same as that in the panel above. The cross 
marks the position of the photometric centre; 



0^2000 = 09^ 26"^ 12.61P, J; 



^02000 



-76° 37' 37.80'^ 



as determined bv lElson et all (|2010[ ). The colour 
scale is described in units of MJy ster~^ by the 
colour bar. 
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Fig. 2.— Top panel: NGC 1705 Hi surface den- 
sity map extracted from the new ATCA Hi synthe- 
sis data. The colour scale is described in units of 
Mq pc~^ by the colour bar. The black ellipse repre- 
sents the B-band R25 radius of 1.2 kpc. The filled star 
marks the central position of the super star cluster, 
NGC 1705-1. The hatched circle in the lower left cor- 
ner of the panel represents the half power beam width 
of the synthesised beam. Bottom panel: 3.6 /xm IRAC 
Spitzer image of the old stellar disk of NGC 1705. The 
black ellipse is the same as that in the panel above. 
The cross marks position of the photometric centre; 
0-2000 = 04^ 541" 13 .50^ ^2000 = -53° 21' 39.80'' 
(|Skrutskie et al.ll2006l ). The colour scale is described 
in units of MJy ster~^ by the colour bar. 



3.2. Ultra-violet and infrared imaging 

Both NGC 2915 and NGC 1705 have been ob- 
serve d with the GALEX satellite ('Martin e t al.l 
2QQ5I) as part of the GALE X Nearby Galaxy Sur- 



) ) as 

vey (|Gil de Paz et al. 20071 ). and with the Spitzer 



satellite as part of the Spitzer Nearb y Galaxies 
Survey f SINGS. iKennicutt et aDl2003[ ). For each 
galaxy, the GALEX far ultra-violet (FUV) imag- 
ing and the Spitzer 24 /im imaging is used to trace 
the directly observable and dust-obscured SFRs, 
respectively. These two tracers are combined to 
yield an estimate for the total SFR (Sec. 14. ip . 

The GALEX FUV band, centred at 1450 A, 
covers the wavelength range 1350 - 1750 A. The 
angular resolution for this band is ^ 4.5''. The 
FUV data were photometrically calibrated and 
corrected fo r attenuation due to dust in the 
Galaxy. The Schlegel et al." (l998) dust extinction 
maps were used to estimate E{B — V) ^ 0."^275 
and E{B - V) - 0."^008 for NGC 2915 and 
NGC 1705, respectively. These B — V extinction 
measures were converted to FUV extinction mag- 
nitudes of Afuv = 8.24 X E{B - V) = 2. ""26 and 
Afuv = 8.24 X E{B - V) = 0.^06 (Wyder et_al] 
20071 ). Neither the FUV nor 24 fim images of 
either galaxy suffer from significant star crowd- 
ing effects. No foreground star subtraction was 
performed. 

The MIPS instrument onboard the Spitzer 
satellite has a field-of-view of 5.4' x 5.4' in the 
24 /im band with a re solution of ^ 6" a nd high 
signal-to-noise ratios. iBidel et al.l (|2008h allude 
to the fact that the MIPS point-spread-function 
at 24 jam is severely non-Gaussian, yet presents no 
significant problems when the data are smoothed 
to resolutions of ~ 20". Again, no foreground star 
subtraction was performed for either galaxy. 

To facilitate a direct comparison of the FUV, 
24 /im and H I data of each galaxy, all images were 
placed on the same astrometric grid with the same 
spatial resolution. The FUV and 24 /im images of 
each galaxy were smoothed, using a Gaussian con- 
volution function, to a resolution equal to the full- 
width-at-half maximum of the synthesised beam 
of the corresponding Hi data. Next, these images 
were re-gridded to the same pixel size and astro- 
metric grid as the H I line data. These smoothed, 
re-gridded FUV and 24 /im images of NGC 2915 
and NGC 1705 are shown in Figs. [3] and HI respec- 



Mq yr ^ kpc ^ 
0.01 0.02 0.03 




30^ 20^ 9^26""! 0^0^ 25^50^ 30^ 20^ 9^26""! 0^0^ 25"'50^ 30^ 20^ 9^26"'1 O^O"" 25"'50^ 
R.A. (2000.0) R.A. (2000.0) R.A. (2000.0) 




28"" 27"" 9'26'' 25"" 24' 
R.A. (2000.0) 



Fig. 3. — FUV, 24 /im and total star formation rate surface density maps for NGC 2915. Top left panel: 
GALEX FUV image of the stellar disk. Top middle panel: Spitzer 24 jam image of the stellar disk. The 
FUV and 24 /im images have been smoothed and re-gridded. Top right panel: star formation rate surface 
density map for the stellar disk, constructed by combining the FUV and 24 /im images according to Eqn. O 
Bottom panel: The full star formation rate surface density map. The single black contour marks the edge 
of the Hi disk. The white contour within the stellar disk is at a level of 0.0018 yr~^ kpc~^. 
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Fig. 4. — FUV, 24 /im and total star formation rate surface density maps for NGC 1705. Left panel: 
GALEX FUV image of the stellar disk. Middle panel: Spitzer 24 fim image of the stellar disk. The FUV 
and 24 /im images have been smoothed and re-gridded. Right panel: star formation rate surface density map 
constructed by combining the FUV and 24/im images according to Eqn. [2l The outer black contour marks 
the edge of the Hi disk. The inner black contour within the stellar disk is at a level of 0.1 yr~^ kpc~^. 



tively. 



surface density, T^sfr, in units of yr ^ kpc ^: 



4. Derived data 

The aim of this work is to understand which 
properties of the neutral ISM regulate the star for- 
mation activity within each galaxy. In order to 
do so we require quantitative measures of the star 
formation activity as well as parameterisations of 
the ISM kinematics. This section describes our 
preferred method of deriving total star formation 
rate surface density estimates for each galaxy, and 
also the rotation curve and ISM velocity dispersion 
parameterisations required for the star formation 
models. 

4.1. The total star formation rate 



The prescription of iLerov et al.l (|2008l ) is used 
to linearly combine the FUV and 24 /im imaging of 
each galaxy to yield an estimate of the total SFR 



Mq yr~^ kpc~^ 



■^24 fim 



^0 yr ^ kpc ^ 

^FUV 



yr 1 kpc 2 



(1) 



= 3.2 X 10~' 



24 



MJy ster-i 

where I24 and Ifuv are the 24 jam and FUV sur- 
face brightnesses, respectively. Very importantly, 
[Eerov et al.l (j2008j) show that this choice of co- 
efficients yields T^sfr values that are consistent 



with t he Ha-24 /im calibration of ICalzetti et al 



(|2007[ ) and that when the 24 /im emission is ig- 
no red the FUV T^sf r estimates reduce to those 
of ISalim et"all (|2007l). The calib r ation of T^sfr 



uses the IMF from Calzetti et al. (120071) which is 



a Kroupa-type two-componen t IMF that extends 
to 120 Mq. Some authors fe.g. lMeurer et al.ll2009i ) 
argue that the IMF is not universal, and that the 
slope of the upper end varies systematically. A 
non-universal IMF implies that the SFR measured 
in a galaxy is highly sensitive to the tracer used 



7 



in the measurement. However, in this work we 
assume a standard two-component Kroupa IMF. 
The ^sFR maps for NGC 2915 and NGC 1705 are 
shown in Fig. [3] and Fig. HI respectively. 

4.2. Rotation curve parameterisations 

The star formation models investigated in this 
work require a measure of the galaxy's rotation 
curve as input. To this end the rotation curves 
are parameterised using the function 



V{r) = Vfiat[l-eM-r/lf 



lat) 



(3) 



where Vfiat and / fiat approximate the asymptotic 
velocity of the outer rotation curve and the length 
scale over which it approaches this constant veloc- 
ity, respectively. This parameterisation is easily 
analytically different ia ble. For NGC 2915 , the ro- 
tation curve derived bv lElson et ^ (|2010h is used. 
This rotation curve was derived by fitting a tilted 
ring model to the Hi velocity field. The best- 
fitting parameters are Vfiat = 83.9 km s~^ and 
I flat = 74.8'' (1.5 kpc). The rotation curve for 
NGC 1705 was generated by parameterising the 
Hi line profiles of an integrated position- velocity 
slice extracted from the new Hi data cube (the 
full details of the derivation of this rotation curve 
will be presented in our forthcoming publication 
Elson et al; 2012, in preparation). This rotation 
curve is very s imilar to that derived previously by 



[m eurer et al.l ([l 



, and is best parameterised by 
Vfiat = 72.8 km s-i and Ifiat = 52.9'' (1.3 kpc). 
The rotation curves together with their respective 
parameterisations are shown in FigO 

4.3. H I velocity dispersions 

The star formation models also require an es- 
timate of the gas velocity dispersion. For some 
of the models considered in this work 



Kennicutt 



1983h . lKennicut t7l989') andlMartin & Kennicutt 



(|200ll ) used a constant gas velocity dispersion of 
o'gas = 6 km s~^, without distinguishing between 
the various thermal phases of the ISM. Several au- 
thors have investigated and emphasised the impor- 
tance of a cold phase of the IS M for the stability 
of disk galaxies. As an example, de Blok fc Walter 
(jjoOG) studied the H i line profiles of NGC 6822 to 
identify separate warm and cold phases with veloc- 
ity dispersions of ^ 8.2 km s~^ and ~ 4.4 km s~^, 
respectively. They showed that the velocity dis- 
persion of the cold component, when used with 
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Fig. 5. — Observed rotation curves of NGC 2915 
and NGC 1705 (black- filled circles) and their pa- 
rameterisations (solid black curves). Error bars in 
the top panel represent the r.m.s. spreads in the 
velocities within thin annuli of the fitted tilted ring 
model. Error bars in the bottom panel represent 
the absolute difference in rotation velocity derived 
separately for each side of the galaxy. 



8 



a Toomre Q criterion, gives an optimal descrip- 
tion of ongoing star formation in NGC 6822, supe- 
rior to that using the more conventional dispersion 
value of 6 km s~^ or the 8.2 km s~ ^ value of the 
warmer component. ISchayd (|2QQ4l ) showed that 
the drop in thermal velocity dispersion associated 
with the transition from the warm (T ^ 10"^ K) 
to the cold (T - 10^ K) phase of the ISM triggers 
gravitational instability on a wide range of scales. 

Throughout this work we assume a constant 
value of a gas = 5.0 km s~^ for the velocity dis- 
persion of the cold gas phase of the neutral ISM. 
Adopting a radially-varying velocity dispersion 
based on the Hi second-order moments would be 
difficult in the cases of NGC 2915 and NGC 1705 
due to the fact that near their centres their Hi 
kinematics are not virialised. 

5. Star formation models 

This section deals with the various star for- 
mation models generated for NGC 2915 and 
NGC 1705. Each model is introduced and its 
method of implementation discussed. The results 
ar e presented t ogether with a comparison to those 
of Leroy et al.l (2008) who carried out similar anal- 
yses for their sample of 23 nearby star-forming 
galaxies from THINGS. This is done so that the 
results can be interpreted in the context of other 
typical star- forming dwarf galaxies. 

5.1. Single- fluid Toomre Criterion 

5.1.1. Introduction 



turbances in regions where the Toomre parameter. 



Safronov (|l960h was the ffist to show that per- 



turbations in a thin, rotating gaseous disk can be- 
come unstable to gravitational collapse due to the 
effects of the s elf-gravity of th e disk. Buildin ; 
on the work of ISafronovl (|l960h . iToomrd ^96. 
showed that a stellar disk that is fairly smooth or 
uniform, and that is rotating in approximate equi- 
librium between its self-gravitational and centrifu- 
gal forces, cannot be entirely stable against the 
tendency to gravitationally collapse. This is the 
case unless the random motions within the disk 
are sufficiently large. Today the Toomre Q crite- 
rion is most commonly used to quantify the grav- 
itational growth of perturbations within a thin, 
rotating gaseous disk. According to this criterion, 
the disk should be unstable to axisymmetric dis- 



ttGE 



(4) 



gas 



is less than unity. The self-gravity, pressure and 
kinematics of the gas disk are represented by S^as, 
a gas and a^, respectively. G is the gravitational 
constant. The epicyclic frequency, z^, is a mea- 
sure of the Coriolis or centrifugal forces stemming 
from the rotation of the perturbations. Following 
Kennicuttl ( 1989[ ). k is estimated as 



n{R) = 1.4l| 



RdV 



(5) 



where V is the rotation velocity at a radius R. It is 
the combination of random motions and centrifu- 
gal accelerations from epicyclic eddies that sup- 
ports the ISM against collapse due to self-gravity. 
The Toomre criterion therefore describes the abil- 
ities of perturbations to rotate around their centre 
of gravity and thus their stability against gravita- 
tional collapse. 

Kennicuttl ( 1989h used a sample of 15 spiral 



galaxies to observationally test the Toomre cri- 
terion on galactic length-s cales. More recently, 
Martin fc Kennicutt (1200 ll ) applied the prescrip- 
tion of iToomrd ([196 47 to a sample of 32 star- 
forming spiral galaxies and compared the radial 
distributions of Qgas to Ha radial surface bright- 
ness profiles. For a sub-sample of 26 galaxies with 
well-defined thresholds they found a median value 
of aq = 0.69 ± 0.2 for the ratio of the gas sur- 
face density to critical surface density defined by 
Eqn. m Scrit = (Jgasf^/^G- Using this ratio as 
an empirical calibration for the Toomre criterion, 
Eqn. m becomes 



Q 



gas 



O^QCTgasI^ 
irGT^aas 



(6) 



This empirical version of the Toomre stability cri- 
terion is used throughout this work. Gravitational 
instability of the gas is predicted by Qgas < 1- 

5.1.2. Methodology 

The Toomre Qgas parameter was calculated for 
every resolution element in the Hi surface den- 
sity map of each galaxy. In assuming a con- 
stant H I velocity dispers ion, we hav e emulated th e 
approaches adopted by iKennicuttI (|l989l . Il998h : 
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Martin fc Kennicuttl (|2QQl[ ): iLerov et al.l (|2QQ8h . 

The epicyclic frequency, hz^ was determined us- 
ing Eqn. [5l The parameterisation of the observed 
rotation curve of each galaxy was used to ana- 
lyticahy determine An azimuthahy averaged 
approximation of the epicyclic frequency was thus 
used for each of the resolution elements at a partic- 
ular galactocentric radius. The Qgas maps shown 
in this section can therefore be thought of as the 
Hi map modulated by a radial function that is 
largely determined by the rotation curve. More 
sophisticated maps can be produced by allowing 
the H I velocity dispersion to vary radially, but in 
this work we assume a constant velocity disper- 
sion. 

5.1.3. Results and discussion 

The results for NGC 2915 and NGC 1705 are 
shown in Figs. [6] and [3 respectively. Most strik- 
ing is the fact that, using the above-mentioned in- 
puts, the Toomre criterion predicts the H I disk of 
both galaxies to be sub-critical, that is Qgas > 1 
for every resolution element in the Hi total in- 
tensity map. This result is inconsistent with the 
star-formation activity observed at the centre of 
each galaxy. The lower panel of Fig. [6] shows the 
logarithmic Qgas distribution for NGC 2915 which 
is approximately Gaussian, peaking at log Qgas ^ 
0.55, far-removed from Qgas = 1 required for grav- 
itational instability. 

Despite there being no Qgas < 1 in the pres- 
ence of clear star formation activity at the cen- 
tre of NGC 2915, the Qgas map does still exhibit 
the general structure that one might expect if it 
were to correctly predict the star formation activ- 
ity. The Qgas radial profile (Fig. [6l middle panel) 
shows the lowest Qgas values to occur at the centre 
of the galaxy, as well as a clear drop in the Qgas 
values at the edge of the stellar disk {R = i?25)- 

The situation is similar for NGC 1705 whose 
logarithmic Qgas distribution is double-peaked 
with Qgas > 1 throughout the Hi disk. This 
galaxy's intense star formation activity is not well- 
described by the Qgas criterion. The very large 
estimates of Qgas ^ 10^'^ for the outer H I disk ar e 
consistent with the results of Me urer et al.l (|l998[ ). 
As above, the lowest Qgas values again occur near 
the galaxy's central stellar disk despite the model 
predicting the H I disk to be sub-critical. 
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Fig. 6. — Top panel: Qgas map for the Hi disk 
of NGC 2915. Middle panel: Radial profile of the 
Qgas map. Bottom panel: Distribution of Qgas 
parameters. The Qgas colour scale in the upper 
panel is descrbied by the colour bar. The gas 
is expected to be unstable to gravitational col- 
lapse in regions where the Qgas< 1- The single 
white star formation rate surface density contour 
in the upper panel is at a level of T^sfr = 0.0018 
M0 yr~^ kpc~^. The dotted vertical line shown in 
the middle panel represents the i?-band R25 radius 
at 1.9 kpc. 
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5.1.4- Comparison to other star-forming dwarfs 
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Fig. 7. — Top panel: Qgas map for the H I disk 
of NGC 1705. Middle panel: Radial profile of the 
Qgas niap. Bottom panel: Distribution of Qgas pa- 
rameters. The Qgas colour scale in the upper panel 
is described by the colour bar. The gas is expected 
to be unstable to gravitational collapse in regions 
where the Qgas< 1- The single black star for- 
mation rate surface density contour in the upper 
panel is at a level of T^sfr = 0.1 yr~^ kpc~^. 
The dotted vertical line shown in the middle panel 
represents the 5-band R25 radius at 1.2 kpc. 
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How do the results of NGC 2915 and NGC 1705 
compare to those of other star-forming dwarfs? 



Lerov et al.l ( 2008 ) found almost no area of the in- 
ner disks of their THINGS galaxies to be formally 
unstable to gravitational collapse. Rather, they 
found that Qgas ~ 4 is typical for the region R < 
O.SR25 and that Qgas ^ 10 is common for larger 
radii. In this sense, the results for NGC 2915 
and NGC 1705 are consistent with the sample of 
iLerov et aP (j2008l ) . The authors find no clear evi- 
dence of a Qgas threshold that can unambiguously 
distinguish between star- formin g regions of high 
and low efficiency. iLerov et al.l (|2008) did, how- 
ever, use a gas = 11 km s~^ as their preferred mea- 
sure of the kinetic support of the ISM. They go on 
to point out that under various assumptions re- 
garding the H2 content and the thermal pressure 
of the galaxies, the median value of Qgas in the 
outer disks of dwarfs agrees quite wel l with the 
threshold v alue of Qnn.s determined by iKennicutt 
(|l989h and iMartin fc KennicuttI (|200lh . Our re- 



sults for NGC 2915 and NGC 1705 as weh as those 



of Lerov et al. 


_(2008|) 


Hunter et al. ( 


19981) ^ 



galaxies, found Qgas to be systematically higher 
by a factor of ~ 2 than the value of determined by 
KennicuttI ^998). 



5.2. Stars+Gas Two-Fluid Toomre Crite- 
rion 

5.2.1. Introduction 

The Toomre criterion discussed in the previous 
section incorporates only the gravitational poten- 
tial of the gaseous disk. In addition to its gaseous 
potential, the stellar potential of a galaxy plays 
an important role in regulating the gravitational 
stability of the ISM by contributing to its self- 
gravity. For this reason two-fiuid instability mod- 
els incorporating the gaseous and stellar potentials 
are particularly relevant to systems in which the 
stellar and gas masses are comparable. NGC 2915 
and NGC 1705 have stellar masses of ~ 3.2 and 
8.8 X 10^ M0, respectively. Each galaxy also con- 
tains at least 10^ M© of H I . This section there- 
fore deals with two-fiuid instability models that 
incorporate the stellar and gaseous potentials of 
NGC 2915 and NGC 1705. 



Rafikovl (|2QQll ) determined the instability con- 
dition for a thin, rotating disk composed of 
gaseous and stehar components to be 



1 



1 



-Rn 



1 + q^Rl 



> 1- (7) 



In this equation = /^ct^^^/ttGI]*, q = ka^^r/^^ 
with k being the wave number of the instabihty. 
The radial velocity dispersion of the stars in the 
plane of the stellar disk is represented by <j*,r, and 

Ra — ^ gas / ^^^r • 

5.2.2. Methodology 

The gravitational effects of the stellar mass 
component of each galaxy were determined using 
Spitzer IRAC 3.6 /im imaging. Ellipses fitted to 
3.6 fim surface brightness isophotes of NGC 2915 
suggest a constant inclination of ~ 55° and a posi- 
tion angle of ~ 306° for the old stellar population. 
A radial profile (Fig. [51 top panel) was derived by 
azimuthally-averaging the stellar surface densities 
in thin annuli. The parameterisation of this profile 
is 



log 



10 



-0.012- 



R 



arcsec 



2.21. (8) 



An identical procedure was carried out to gen- 
erate a stellar surface density radial profile for 
NGC 1705 (Fig. El bottom panel). A constant in- 
clination and position angle of 45° was assumed. 
The parameterisation of the resulting profile is 



log 



10 



Mqpc- 



= -0.015- 



R 



arcsec 



2.37. (9) 



All stellar surface densities were inclination cor- 
rected so that Eqns.[8]and[9]parameterise the face- 
on values. 

The two-fluid stability criterion of iRafikov 



(|200ll ) requires an estimate of the radial com- 
ponent of the stellar velocity disper sion of each 
galaxy, a^^r- The prescription of iLeroy et al. 
(|2008l ) was followed to deter mine (j^jfj^T a s a func- 



tion of galactocentric radius. iLeroy et al.l (|2| 



made four assumptions in order to determine the 
radial stellar velocity dispersion: 

1. The scale height of an exponential stellar 
disk is constant with galactocentric radius. 
This is a typical observat ion for edge-on 
galaxies (|Kregel et al.ll2002[ ). 




100 



200 




Fig. 8. — Stellar surface density radial profiles 
derived from the Spitzer IRAC 3.6 jiia imaging 
of each galaxy (black filled circles). All surface 
densities have been inclination corrected. The pa- 
rameterisation of each profile is shown as a solid 
black line. The dashed vertical lines mark the i?25 
isophotal radii of the galaxies. 
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2. A flattening ratio of U/h^ = 7.3 ± 2.2, 
where and are the disk scale length 
and scale height, respectively. This is the 
average volume-corrected flattening of the 
disk light for the s ample of 34 edge-on spi- 
ral galaxi es of | de Gr iis (1998) as determined 
by iKregel et all (|2QQ2f ). From the IRAC 
3.6 fim imaging the stellar disk scale lengths 
of NGC 2915 and NGC 1705 were deter- 
mined to be U ~ 0.6 kpc and U ~ 0.14 kpc, 
respectively, leading to ~ 0.08 kpc and 
K ^ 0.02 kpc. 

3. The stellar disk is isothermal in the z- 
direction. Under the assumption of hy- 
drostatic equilibrium, this leads to /i* = 

1/2 J cri ^/ 2110 (|van der Kruit fc Searle 



Il98ll ). where is the mid-plane stellar vol- 
ume density. Using = Ap^h^ and elim- 
inatin^ then yi elds a^^z = V^ttGS^/i* 
( van der Kruit|[l988 ). leading to 



27iGl^Yi^ 



7.3 



(10) 



4. A flxed ratio of cr^^^/c^*,^ 



0.6 to relate 



the vertical and radial stellar velocity dis- 
persions. This assumption is based on lim- 
ited available eviden ce for late-type systems 
(IShapiro et all [2 0031 ). 



The parameterisation of the stellar surface den- 
sity radial proflle of each galaxy was used together 
with the approximation to yield cr^^r as a func- 
tion of galactocentric radius. Figure [9] displays 
these a^^r radial proflles. At inner radii they are 
similar to the constan t a^^r = 25 km s~^ estimate 
used bv lRaflkovl ( 200ll ) for the implementation and 
testing of his two-fluid stability criterion. The ra- 
dial stellar velocity dispersion monotonically ap- 
proaches zero beyond the stellar disk. 
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Fig. 9. — Radial proflles of the stellar velocity dis- 
persion parallel to the plane of the stellar disk 
for NGC 2915 and NGC 1705 (solid and dashed 
curves, respectively), obtained by parameterising 
the stellar surface density radial proflles. 
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So far, the prescriptions used to estimate 
and cr*,r have been presented. These parameters 
together with the epicychc frequency radial pro- 
file of each galaxy allow Q* in Eqn. [Tito be deter- 
mined. The other parameters in Eqn. [3 that need 
to be considered are q and Ra- The q parame- 
ter is a function of the wave number k = 27r/A, 
of the inst ability. Specifying k determines q. Ac- 
cording to iRosolowskv fc Blitzl (j2QQ5) . 80% of the 
Milky Way's molecular gas is contained within gi- 
ant molecular clouds ranging in size from ~ 2 pc 
to ~ 100 pc. Since we have no a priori knowl- 
edge of the linear size of a typical giant molecular 
cloud in NGC 2915 or NGC 1705, a set of pertur- 
bation wave numbers was selected that spans this 
length scale range. Perturbation length scales of 
A = 0.005, 0.1, 0.2, 0.3, 0.5, 0.8 kpc were con- 
sidered, each resulting in a separate instability 
map. These perturbation length scales correspond 
to wave numbers of k = 40, 20, 10, 6, 4, and 2.5 
times TT. In the context of the instability maps, 
the length scales we consider range in linear extent 
from that of a single resolution element to that of 
the entire gaseous disk. Finally, Ra = (Jgas/o'^^r 
was specified by assuming a gas = 5.0 km s~^ for 
all two-fiuid instability maps and using the cr*,^ 
profiles shown in Fig. [9l 

5.2.3. Results and discussion 

The results of the 1/Q^as,* two-fiuid instability 
analyses for NGC 2915 and NGC 1705 are pre- 
sented in Figs. [To] and [TTJ The radial profiles of 
these maps are shown in Fig. [12] while the distri- 
bution of 1/Qgas,* is shown in Fig. [T3l 

The NGC 2915 results are considered first. 
Even with the gravitational potential of the stel- 
lar disk included, the galaxy is sub-critical with 
< 1 throughout its Hi disk for all 
implemented wave numbers. This two fiuid in- 
stability criterion, like the single fiuid instability 
criterion considered in the previous section, for- 
mally fails to predict the observed star formation 
activity at the centre of NGC 2915. The extra 
self-gravity of the ISM induced by the stellar po- 
tential is still not enough to counteract the large 
epicyclic frequencies at inner radii. Despite this 
formal failure for the central part of the galaxy, the 
instability maps again produce the sort of struc- 
ture that one might expect for a late-type galaxy. 
Beyond the radial extent of the stellar core, the 



two-fiuid l/Qgas,* maps resemble the single-fiuid 
Qgas maps. This result is expected since the stel- 
lar potential is negligible in the outer gaseous disk 
of the galaxy. Including the gravitational infiu- 
ence of the stellar disk of NGC 2915 leaves the 
inner portion of its HI disk only marginally sta- 
ble against gravitational collapse for the small- 
est wave numbers, corresponding to perturbation 
length-scales of ~ 0.5 kpc. This r e sult is fairly 
consistent with that of iLerov et al. ( 2008[ ). The 
authors found no clear evidence of a single-fiuid 
Toomre Qgas criterion that could unambiguously 
distinguish between high and low star formation 
efiiciency regions. However, when using Eqn. [7] 
to include the effects of stellar gravity they found 
some regions became gravitationally unstable. In- 
cluding the stellar gravity led to large portions 
of the gaseous disks in their sample being only 
marginally stable against large-scale gravitational 
collapse, similar to our NGC 2915 result. 

The situation is only slightly better for the 
NGC 1705 instability maps which for some of the 
implemented wave numbers have 1/Qgas,* > 1- 
The Qgas,* = 1 level is marked with a single black 
contour in the instability maps shown in Fig. [TTJ 
The best instability model is that for a wave num- 
ber of k = IOtt, which corresponds to a pertur- 
bation length-scale of A = 200 pc. This model 
predicts the very central regions of the H I disk to 
be unstable. In the context of a two-fiuid Toomre 
criterion, the gravitational potential of the stel- 
lar disk seems to contribute enough self-gravity to 
the ISM to yield the central disk unstable. Large 
perturbation length-scales of A ~ 0.8 kpc make 
the entire Hi disk of the galaxy unstable against 
gravitational collapse. This demonstrates the in- 
appropriateness of the perturbation size. Pertur- 
bations with wave number k ^ lOir yield only the 
inner disk unstable while the degree of instability 
is roughly constant over the outer Hi disk. The 
f/Qgas,* parameters for NGC 1705 general l y span 
a small range of values. As iLerov et"al] ( 20081 ) 



point out, a small spread in 1/Q^as,* values is con- 
sistent with a self-regulating star formation sce- 
nario in which the gravitational potential asso- 
ciated with newly formed stars will seed further 
gravitational collapse of the ISM. 

Finally, it should be mentioned that several 
sources of uncertainty affect the results of these 
two-fiuid instability analyses. The main source 
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Fig. 10. — NGC 2915 stars+gas two-fluid instability maps for various perturbation length scales. The 
common colour scale is specified by the colour bar at the top of the figure, and represents the two-fiuid 
instability parameter, l/Qgas,*- The perturbation wave number, k = 27r/A, for each instability map is 
shown in the top left corner of the panel. Wave numbers of k = 40, 20, 10, 6, 4, 2.5 times tt correspond 
to wave lengths of A = 0.005, 0.1, 0.2, 0.3, 0.5, 0.8 kpc, respectively. The gaseous disk is expected to be 
unstable to large-scale gravitational collapse in regions where 1/Qgas,* > 1- For all of the above-presented 
maps, the instability criterion predicts a sub-critical gaseous disk for NGC 2915. 
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Fig. 11. — NGC 1705 stars+gas two-fluid instability maps for various perturbation length scales. The 
intensity scale of each map is specified by the colour bar at the top of the panel, and represents the two-fiuid 
instability parameter, l/Q^as,*- The perturbation wave number, k = 27r/A, for each instability map is shown 
in the top left corner of the panel. Wave numbers of k = 40, 20, 10, 6, 4, 2.5 times tt correspond to wave 
lengths of A = 0.005, 0.1, 0.2, 0.3, 0.5, 0.8 kpc, respectively. The gaseous disk is expected to be unstable to 
large-scale gravitational collapse in regions where 1/Q^as,* > 1- The single black contour shown in some of 
the panels is at a level of 1/Q^as,* = I, enclosing the region of the galaxy in which gravitational instability 
is formally expected. 
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Fig. 12.— Radial profiles of the NGC 2915 and 
NGC 1705 stars+gas two-fluid instability maps 
shown in Figs. [10] and [Til respectively. From top 
to bottom the profiles correspond to wave numbers 
oik = 40, 20, 10, 6, 4, 2.5 times tt. 



is the lack of direct measurements of the stellar 
velocity dispersions. The observed stellar surface 
density radial profile of each galaxy has been used 
to infer a stellar velocity dispersion profile. This 
appro ach assumes that a* oc (iLerov et al 



20081 ) . This assumption may well affect the results 



of the star/gas two-fluid stability analyses. 

5.2.4- Comparison to other dwarfs 
A main finding of iLerov et al.l 



was that 

considering the contribution of the stellar poten- 
tial to the self-gravity of the ISM does not ren- 
der large, extended portions of the disk unsta- 
ble. The NGC 1705 re sults are well aligned with 
the iLerov et al.l (|2QQ8h findings, with it having 
a significant fraction of its inner disk being for- 
mally unstable. These results should be com- 
pared to the corresponding results for the sim- 
ple Toomre Q case in which only the gaseous 
potential was considered and for which no for- 
mal gravitational instability was predicted for ei- 
ther galaxy (i.e. Qgas > 1 throughou t the d isk). 



Another main finding by iLerov etliD ( 2QQ8I ) was 
that the 1/Qgas,* parameters exhibit a narrower 
range of values than the Qgas parameters. This is 
clearly the case for NGC 1705 for which 1/Qgas,* 
spans a range of y 0.3 dex for all wave numbers. 
iLerov etliD ( 2008 ) suggest that a narrow range of 



values means that 1/Qgas,* offers little 
leverage to predict the locations of efficient star 
formation. 



5.3. Dark matter+gas two- fluid Toomre 
criterion 

5.3.1. Introduction 

The results of the previous section demonstrate 
the important role that the stellar gravitational 
potential of a galaxy can play in regulating the 
star formation activity. This section attempts to 
answer the question of whether the dark matter 
(DM) can play a similar role. 

NGC 2915 and NGC 1705 are both DM- 
dominated galaxies. NGC 2915 contains about 
1.4 X 10^^ Mq of dark mass within ~ 10 kpc of 
its centre and has its DM strongly concentrated 
with a core density of .17 M^ P^~^ ^ core 
radius of - 0.9 kpc (|Elson et all 12010! ). This 
is for the case in which the halo is modelled as a 
pseudo-isothermal sphere. The situation is similar 
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for NGC 1705 which contains - 3.1 x 10^ of 
DM within ~ 6 kpc of its centre, and which has 
a core density of 0.01 Mq pc~^ and a core ra- 
dius of 1.2 kpc (Elson et al. 2012, in prep.). Since 
these galaxies contain dense concentrations of DM 
where the active star formation is observed, it is 
feasible that the portion of the DM halo of each 
galaxy that is co- located with its H I disk plays a 
role in regulating the star formation activity. 

For NGC 2915, iBureau et all (|l999h proposed 
the self-gravitational effects of a large dark matter 
component co-located with the system's Hi disk 
as a possible e xplanation for its o bserv ed H I spi- 
ral structure. 'Masset fc Bu reau! ( 20031 ) used hy- 
drodynamical simulations to further explore this 
proposed mechanism, comparing the results to the 
observations using customised column density con- 
straints. They found that when the observed Hi 
density is scaled up by a factor of 10, the disk de- 
velops a spiral structure that closely resembles the 
observed one. They suggest the disk of NGC 2915 
to contain a large dark mass component that is 
closely linked to the observable neutral ISM. They 
show further that the proposed scahng does not re- 
sult in wide-spread star formation in the outer H I 
disk, consistent with observations. 



iHunter et al.l (^199&) state that disk dark matter 
"effectively acts like stars in a two-fluid instability, 
giving extra self-gravity to small perturbations in 
the gas" . In this section we therefore test an analo- 
gous version of the stars+gas instability criterion 
from Sec. 15.21 by replacing the stellar quantities 
with corresponding DM quantities. The two-fluid 
DM+gas instability criterion is 

1 2^2 



-Ra 



Qgas.DM QdM 1 + Q gas 1 + Q^R^ 

(11) 

where Qdm = i^ctdm /ttG^dm, q = kaoM/i^, 
with k being the wave number of the instability. 
The velocity dispersion of the DM is represented 
by (TDM, and Ra = (Jgas/crDM' This model, as 
compared to a simple Toomre Q instability model 
with the Hi surface densities boosted to account 
for the disk dark matter, has the advantage that it 
incorporates a separate measure of the dark mat- 
ter distribution and velocity dispersion. The ve- 
locity dispersion of the dark matter is expected to 
be much larger than that of the H I , meaning that 
equal masses of gaseous and dark matter will have 
different effects on the gravitational instability of 



the gas. 

5.3.2. Methodology 

S DM radial profiles for the galaxies are required 
by Eqn. [TTJ Each galaxy has its dynamics best 
explained by a pseudo-isothermal sphere parame- 
terisation of its DM halo: 



pDM{r) 



1 



(12) 



where po is the DM core density and Tc the core 
radius. The DM surface density (in units of 
M0 pc~^) along a particular line of sight is mod- 
elled as the DM volume density at the correspond- 
ing galactocentric radius integrated over the thick- 
ness of the H I disk: 



^DM(r) 



hgas 



PDM{r)dh 



(13) 



po 



hgas 1 + {r/rc) 

2 hgas po 



jdh (14) 



1 + (r/rc) 



2 • 



(15) 



Hi scale heights of hgas = 0.41 and 0.24 kpc are 
adopted for NGC 2915 and NGC 1705, respec- 
tively. These estimates were ge nerated using the 
prescription of Walter fc Brinksl ([l999). Adopting 
the pQ and Vc pseudo-isothermal sphere parame- 
ters mentioned in Sec. l5.3| the T^dm radial profiles 
shown in Fig. [14] are obtained. 

Also required by Eqn. [TT] is an estimate of 
the DM velocity dispersion. In lieu of any such 
determination for either galaxy, we use the 3- 
dimensional velocity dispersion of an isotropic 
isothermal sphere: 



(16) 



(^DM = \l -Vflati 



where Vfiat is the asymptotic velocity of the outer 
rotation curve. Adopting the best-fitting Vfiat pa- 
rameters from Sec. 14. 2[ a dm esitmates of 102.7 
and 89.2 km s"^ are obtained for NGC 2915 and 
NGC 1705, respectively. Finally, whereas a range 
of wave numbers were tested in Sec. 15. 2[ here we 
test single wave numbers of k = 2.57r and 47r for 
NGC 2915 and NGC 1705, respectively. 
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5.3.3. Results and discussion 




-0.5 

logio(Qgas,*) 

Fig. 13.— Distributions of the NGC 2915 and 
NGC 1705 stars+gas two-fluid instability maps 
shown in Figs. [TQland [TT| respectively. 




100 200 300 
R (arcsec) 

Fig. 14.— NGC 2915 and NGC 1705 dark matter 
surface density profiles (solid and dashed curves, 
respectively), as specified by Eqn.[T5l used to con- 
struct the gas+DM two-fluid instability maps. 



The results for NGC 2915 and NGC 1705 are 
presented in Figs. [15] and [161 respectively. The 
results suggest that part of the inner Hi disk of 
each galaxy can become unstable to gravitational 
collapse when additional self-gravity of the gas 
caused by the disk dark matter is incorporated. 
Disk dark matter could therefore contribute to- 
wards regulating the star formation activity near 
the centres of NGC 2915 and NGC 1705. 

Using each galaxy's pseudo-isothermal sphere 
parameterisation to estimate the amount of DM 
co-located with its Hi disk ensures that the in- 
ner Hi disk (where pDuij) is high) contains more 
dark mass than the outer Hi disk (where pDuij) 
is relatively lower). The resulting distribution of 
^ DM values is different to the one obtained by sim- 
ply boosting each line-of-sight T^hi measurement 
by a constant factor. 

I n their review of col d gas accretion in galax- 
(|2008) 



les 



Sancisi et al 



allude to the fact that 
"it is remarkable that there is such a pronounced 
spiral structure in the outer regions of spirals 
where dark matter dominates and even in dwarfs 
where the dark halo is believed to be predomi- 
nant everywhere". The question, they therefore 
say, is "whether these systems have light disks 
surrounded by massive dark halos or, rather, have 
heavy and dark disks" . The modeling results from 
this section allow this question to be partly ad- 
dressed for the cases of NGC 2915 and NGC 1705. 
It has been shown that only a heavy disk sce- 
nario can account for the observed star formation 
in both NGC 2915 and NGC 1705. This result, to- 



gether with the result of Masset fc Bureaul (jli 



that a heavy disk allows for the natural formation 
of NGC 2915's observed Hi spiral structure sug- 
gests that certainly NGC 2915 may have a heavy 
Hi disk, and very likely NGC 1705 too, with sig- 
nificant amounts of dark mass distributed within 
each of them. 

5.4. Shear-based instability criterion 

A common feature of all the models pre- 
sented so far is their characterisation of the global 
kinematics of a galaxy-they all incorporate the 
epicyclic frequency. This section investigates a 
star formation model for which the kinematics are 
quantified not in terms of the Coriolis force (i.e. 
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Hi), but rather the rotational shear. 

5.4- Introduction 

iHunter et aP (|l998h suggest the Coriohs force 
incorporated in the Toomre Q instabihty criterion 
to become less important in the presence of rota- 
tional shear. It is the time available for clouds to 
collapse in the presence of rotational shear, they 
argue, that regulates a gal axy's star formation ac- 



tivity. iHunter et al.l () 19981 ) use Oort's A constant, 



r ar J 



(17) 



to quantify the rotational shear of the gas. Their 
so-called shear-based parameter for gravitational 
instability is then 



aA(TgasA 

^gas — 

TTUljgas 



(18) 



Based upon the idea that a perturbation must 
grow by a factor of ~ 100 in the pre sence of shear 



for th e instability to be significant, IHunter et al 
(|l998l) estimate a a ~ 2.5. This value of A matches 
the contrast between the surface densities of neu- 
tral and molecular inter-stellar media in the pres- 
ence of rotational shear. It also allows the condi- 
tion Qgas ^ 1 to be met over the V oc R portion 
of the rotation curve. Regions within the galaxy 
that have Sgas < 1 should be forming stars while 
regions with Sgas > 1 should be stable against 
large-scale gravitational collapse. 

'Hunter e t al.l (|l99 8l) point out that using A to 
quantify the gas kinematics instead of makes 
very little difference for a flat rotation curve. They 
estimate the two thresholds to be the same to 
within 12% for a flat rotation curve. The differ- 
ence between the two thresholds is most appar- 
ent for rising rotation curves where the rotational 
shear is low. Under such circumstances Sgas can 
be significantly lower than Qgas • Figure [T71 shows 
radial profiles of the epicyclic frequency and rota- 
tional shear for the parameterisation of the rota- 
tion curve of NGC 2915 (Fig.[5j). Most noticeable 
is the fact that the rotational shear is a factor 
~ 5 — 30 smaller in magnitude than the epicyclic 
frequency within R < 150'^ The situation for 
NGC 1705 is very similar, and it is therefore within 
the spatial extent of the stellar disks of the galax- 
ies that a the shear-based characterisation of the 



kinematics is expected to be most influential when 
incorporated into a star formation model. 

5.4.2. Methodology 

Sgas maps were produced for NGC 2915 and 
NGC 1705 according to Eqn. [HI As for in 
the Toomre Qgas criterion, the shear parameter, 
A, for each line of sight was approximated as the 
azimuthally- averaged shear corresponding to the 
galactocentric radius of the resolution element. 
The shear radial profile was determined for each 
galaxy using Eqn. [iTl together with the parame- 
terisation of the galaxy's rotation curve. 

5.4-3. Results and discussion 

The Sgas maps for NGC 2915 and NGC 1705 
are shown in Figs. [18] and [191 respectively. The 
Sgas criterion correctly describes the star forma- 
tion activity at the centre of each galaxy's gaseous 
disk. The Sgas = 1 contours of the shear maps are 
able to accurately trace the edges of the stellar 
disks of the galaxies. By simply using a shear- 
based characterisation of the global dynamics, the 
Sgas criterion can correctly locate the unstable 
portions of the disks without the need for extra 
stellar or dark matter self-gravity. There appears 
to be no tight correlation between Sgas and the ob- 
served star formation activity. Within the central 
star-forming region of each galaxy, the Sgas pa- 
rameters vary significantly (e.g. a factor of ~ 10 
for the case of NGC 1705's stellar disk). 

5.4-4- Comparison to other dwarfs 

The NGC 2915 and NG C 1705 f indings are con- 
sistent with the results of iLerov et al. (2008) who 
found the inner disks of their galaxies to be more 
nearly super-critical in the context of a shear- 
based star formation threshold than in the con- 
text of a Qgas criterion. Consistent with our re- 
sults foLiiG^915_aiid NGC 1705 is the fact 
that Leroy et al.l (|200 8l) find a simple shear cri- 
terion to perform better than the star+gas two- 
fluid criteri^n_for their sample of THINGS galax- 
ies. Leroy et al.l (|2008[ ) found that, according to 



this shear-based instability criterion, there are re- 
gions of high star formation activity within the 
disk that are formally unstable against the ten- 
dency to gravitationally collapse, instead of being 
formally stable. 
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6. Molecular hydrogen content 

Nowhere in this work have we made provi- 
sions in the star formation models for a possible 
molecular component of each galaxy's ISM. Stars 
are known to form from molecular clouds, and 
so the presence of molecular gas is expected for 
each galaxy Despite this expectation, however, 
no published CO detections exist for NGC 2915 
or NGC 1705. In this work we therefore could 
not directly study the possible links between the 
star-forming properties of the galaxies and their 
molecular inter-stellar media. 



Bigiel et al. ( 2QQ8[ ) showed that a Schmidt-type 
power law with index N = 1.0 ± 2 relates T^sfr to 
the molecular gas surface density {^H2) across a 
sample of THINGS spiral galaxies, implying that 
H2 forms stars at a constant efficiency. This re- 
lation allows an estimate of the molecular con- 
tent of a system to be obtained from its observed 
star formation activity. Using the total star for- 
mation rate surface density maps of NGC 2915 
and NGC 1705, respective total molecular gas 
masses of Mh^ = ^-^-I'.i ^ ^0 ^nd = 



7.4I7 Q X 10^ Mq are inferred to be contained 



within each galaxy's R25 isophotal radius. The H I 
masses of NGC 2915 and NGC 1705 are k nown 
to be - 5.5 X 10^ M.:. (lElson et al.l l2010h and 
1.1 X 10^ M0 (jMeurer et al.lll998h . respectivelv. 
suggesting molecular-to-atomic-H I mass ratios of 
Mh2/Mhi 0.08 and Mh2/Mhi 0.48. These 
mass ratios are lower than the average ratio of 
Mh2/Mhi = 0.89 for the 23 gas-ri ch late-type 



galaxies from the THINGS sample of iLerov et al 



( 20081 ). This is not surprising given the small sizes 
of the stellar disks of NGC 2915 and NGC 1705 
relative to the sizes of their extended H I disks. 

Radial profiles of the SFR surface density maps 
were used together with the above-mentioned 
Schmidt-type power law from [Bigiel et ah! (2008) 
to generate the H2 surface density radial pro- 
files shown in Fig. [20l For comparative purposes 
each galaxy's Hi profile is also shown. The in- 
ferred H2 surface densities within the R23 ra- 
dius of NGC 2915's stellar disk are much lower 
(< 2 M0 pc~^) than the corresponding Hi sur- 
face densities. Incorporating them into the star 
formation laws would result in an increase of the 
gas surface density by a factor of 1.2 at most. Such 
an increase will not significantly affect the results 



presented in this work. From the Qgas radial pro- 
ffies of NGC 2915 and NGC 1705 presented in 
Figs. [6] and [71 it is clear that a boost factor of at 
least ~ 2 is required to yield the inner most por- 
tions of each galaxy's gaseous disk gravitationally 
unstable. 

Because of the intense star burst activity near 
the centre of NGC 1705, some of the system's in- 
ferred H2 surface densities are much higher than 
the corresponding Hi surface densities. Caution 
has to be exercised in the case of NGC 1705, 
however. Altho ugh the THINGS s ample of late- 
type galaxies of iBigiel et al.l ( 2008h does include 
11 late- type dwarfs, none of those systems ex- 
hibit star formation activity as extreme as that 
of NGC 1705. It is highly likely that extreme sys- 
tems such as NGC 1705 do not obey the same 
star formation law as do more typical star form- 
ing systems. In this sense we cannot quantitatively 
comment on the extent to which an incorporation 
of NGC 1705's H2 content into the star forma- 
tion models will affect the results presented in this 
work. 

7. Summary 

This paper deals with the investigations of var- 
ious star formation models for the blue compact 
dwarf galaxies NGC 2915 and NGC 1705. Both 
galaxies contain small stellar disks embedded in 
much larger H I disks. We have used GALEX far 
ultra-violet and Spitzer 24 /im imaging of each 
galaxy to estimate its total star formation rate. 
New deep, high-resolution Hi synthesis observa- 
tions of each system have been used to quantify 
the distribution and kinematics of its neutral ISM. 

The star formation models incorporate vari- 
ous properties of each galaxy's ISM to determine 
which regions within the galaxy should be unsta- 
ble against the tendency to gravitationally collapse 
and hence form high-mass stars on varying galac- 
tic length scales. The success or failure of these 
models sheds light on the interplay between the 
key properties of the ISM that regulate the star 
formation activity. The star formation thresholds 
in NGC 2915 and NGC 1705 are not purely local 
phenomena. The Toomre single-fluid Qgas crite- 
rion as well as the stars+gas and DM+gas two- 
fluid criteria {Qgas,* and Qgas, dm respectively) all 
use the epicyclic frequency at a given galactocen- 
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trie radius within the galaxy as a measure of its 
kinematics. The Qgas criterion incorrectly pre- 
dicts the entire Hi disks of both NGC 2915 and 
NGC 1705 to be sub-critical. This result is con- 



sistent with that of iLerov et al.l (|2QQ8l ) who found 
almost no area of the inner disk for their sample 
of 23 THINGS galaxies to be formally unstable to 
gravitational collapse in a Toomre Qgas context. 
These results suggest the Toomre Qgas criterion 
to be of limited utility in terms of distinguishing 
star-forming regions of high and low efficiency. 

The stellar potential of each galaxy is expected 
to play a crucial role in regulating the star for- 
mation activity. Despite this expectation, how- 
ever, the Qgas,* criterion also fails for NGC 2915 
whose entire H I disk is predicted to be sub-critical. 
The model cannot describe NGC 2915's central 
star formation activity. The situation improves for 
NGC 1705 for which perturbation length scales of 
A ^ 200 pc lead to predictions of its inner H I disk 
being formally gravitationally unstable. 

Disk dark matter can account for the ob- 
served star formation activity in NGC 2915 and 
NGC 1705. Both galaxies are known to be ex- 
tremely dark-matter-dominated. The fraction of 
their dark matter that is co- located within the H I 
disk is treated as contributing to the self-gravity of 
the ISM. Assuming a pseudo isothermal sphere pa- 
rameterisation of each galaxy's dark matter halo, 
both NGC 2915 and NGC 1705 have formally un- 
stable inner Hi disks. This result does not hold 
for the case in which the distribution of the as- 
sumed disk dark matter exactly traces the Hi 
distribution. 

The final star formation model investigated in 
this paper is a shear-based model built on the 
premise that it is the time available for clouds to 
collapse in the presence of rotational shear that 
regulates a galaxy's star formation. Indeed, adopt- 
ing the rotational shear as a characterisation of 
the global kinematics allows the Sgas criterion 
to accurately locate the unstable parts of both 
NGC 2915's and NGC 1705's H i disks without the 
need for stellar or dark matter self-g ravity. The 
same general result was obtained by iLerov et al. 



(|2008f ) who found the inner disks of their galaxies 
to be more nearly super-critical in the context of 
a shear-based Sgas criterion than in the context of 
a Qgas criterion. 

In conclusion, despite NGC 2915 and NGC 1705 



both being unusual star-forming galaxies, their 
star-formation activity can be partly understood 
in terms of the self-gravity of the ISM which, in 
turn, is controlled by the combined effects of the 
gravitational potentials of their various mass com- 
ponents. Alternatively, the star formation activity 
in each system can be accounted for by correctly 
characterising its global kinematics. Throughout 
this work, in lieu of published CO data, we have 
ignored any possible H2 components. Molecu- 
lar hydrogen must be present, it is unlikely that 
the observed star formation could be taking place 
without it. Based on our measured star formation 
rates, we have estimated the amount of H2 within 
each galaxy to be of the order of ~ 10 M© {r^ 2—3 
orders of magnitude less than the dynamical mass 
of each galaxy). Although we cannot reliably in- 
corporate the effects of this mass component into 
our star formation models, they are not expected 
to significantly affect the results. 
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Fig. 15. — Gas+DM two-fluid instability map 
for NGC 2915. The colour bar describes the 
^/Q gas, DM valucs. The pseudo-isothermal pa- 
rameters used to generate the map are specified 
in the lower left corner, po and Tc are in units 
of Mq pc~^ and kpc, respectively. The edge of 
the stellar disk is delimited by the solid white 
^SFR = 0.0018 Mq pc~^ contour while the grav- 
itationally unstable portion of the gas disk is out- 
lined by a solid black l/Qgas,DM = 1 contour. 
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Fig. 16. — Gas+DM two-fiuid instability map 
for NGC 1705. The colour bar describes the 
1/ Q gas DM values. The pseudo-isothermal pa- 
rameters used to generate the map are specified 
in the lower left corner, po and Tc are in units 
of Mq pc~^ and kpc, respectively. The edge of 
the stellar disk is delimited by the solid white 
^SFR = 0.1 M0 pc~^ contour while the gravita- 
tionally unstable portion of the gas disk is outlined 
by a solid black l/Qgas,DM = 1 contour. 
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Fig. 18. — Top panel: Sgas instability map for 
the Hi disk of NGC 2915. Middle panel: Ra- 
dial profile of the Sgas instability map. Bottom 
panel: Distribution of logiQ{Sgas) values. The 
colour bar above the upper panel describes the 
Sgas values. The single yellow contour in the map 
is at a level of Sgas = 1, and encloses the unstable 
region of the Hi disk. The single white contour 
in the upper panel, at a level of T^sfr = 0.0018 
M0 yr~^ kpc~^, approximates the edge of the stel- 
lar disk of NGC 2915. 
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Fig. 20. — Observed Hi and inferred H2 radial 
profiles (blue-dashed and solid red curves, respec- 
tively) for NGC 2915 and NGC 1705 (top and bot- 
tom panels, respectively). Note the logarithmic 
ordinate axis of the bottom panel. The H I radial 
profi les for NCG 2915 and NGC 1705 are taken 
from lElson et al.l (|2010[ ) and Elson et al. (2012, in 
prep.), respectively. 



Fig. 19. — Top row: Sgas instability map for the 
Hi disk of NGC 1705. Middle row: Radial profile 
of the logiQ{Sgas) values. Bottom row: Distribu- 
tion of logiQ{Sgas) values. The colour bar above 
the upper panel describes the Sgas values. The 
single black contour in the map is at a level of 
Sgas = 1, and encloses the unstable region of the 
Hi disk. The single white contour in the upper 
panel, at a level of T^sfr = 0.1 Mq yr~^ kpc~^, 
approximates the edge of the star-forming core of 
NGC 1705. 



